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We analyze non - Fermi liquid (NFL) behavior of fluctuating gap model (FGM) of pseudogap 
behavior in both ID and 2D. We discuss in detail quasiparticle renormalization (Z - factor), demon- 
strating a kind of "marginal" Fermi liquid or Luttinger liquid behavior and topological stability of 
the "bare" Fermi surface (Luttinger theorem). In 2D case we discuss effective picture of Fermi 
surface "destruction" both in "hot spots" model of dielectric (AFM, CDW) pseudogap fluctuations, 
as well as for qualitatively different case of superconducting d - wave fluctuations, reflecting NFL 
spectral density behavior and similar to that observed in ARPES experiments on copper oxides. 



PACS numbers: 71.10.Hf, 71.27. +a, 74.72.-h 



I. INTRODUCTION 

Pseudogap formation in the electronic spectrum of un- 
derdoped copper oxides is especially striking anomaly of 
the normal state of high temperature superconductors^. 
Discussions on the nature of the pseudogap state continue 
within two main "scenarios" - that of superconducting 
fluctuations, leading to Cooper pair formation above T c , 
or that of another order parameter fluctuations, in fact 
competing with superconductivity. 

We believe that the preferable "scenario" for pseudo- 
gap formation is most likely based on the model of strong 
scattering of the charge carriers by short-ranged antifer- 
romagnetic (AFM, SDW) spin fluctuations^. In momen- 
tum representation this scattering transfers momenta of 
the order of Q = (-,-) (a — lattice constant of two di- 
mensional lattice) . This leads to the formation of struc- 
tures in the one-particle spectrum, which are precursors 
of the changes in the spectra due to long-range AFM 
order (period doubling). 

Within this spin-fluctuation scenario a simplified 
model of the pseudogap state was studied^S^ under the 
assumption that the scattering by dynamic spin fluctu- 
ations can be reduced for high enough temperatures to 
a static Gaussian random field (quenched disorder) of 
pseudogap fluctuations. These fluctuations are defined 
by a characteristic scattering vector from the vicinity of 
Q, with a width determined by the inverse correlation 
length of short-range order re — Actually, a simi- 

lar model (formalism) can be applied also to the case of 
pseudogaps of superconducting nature^. 

These models originated from earlier one - dimensional 
model of pseudogap behavios^, the so called fluctuat- 
ing gap model (FGM), which is exactly solvable in the 
asymptotic limit of large correlation lengths of pseudogap 
fluctuations re = — > 0», and "nearly exactly" solvable 
case of finite re, where we can take into account all Feyn- 
man diagrams of perturbation series, though using an 
approximate Ansatz for higher - order contributions 5 . 

Non - Fermi liquid behavior of FGM model was dis- 
cussed already in the case of lD^Si^, as well as in 
2£ji i 2 i 3_ However, some interesting aspects of this model 



are still under discussion 9 . Below we shall analyze differ- 
ent aspects of this anomalous behavior both in ID and 
2D versions, mainly for the case of AFM (SDW) or CDW 
pseudgap fluctuations, and also, more brielfly for the case 
of superconducting fluctuations, demonstrating a kind of 
"marginal" Fermi liquid behavior and qualitative picture 
of Fermi surface "destruction" and formation of "Fermi 
arcs" in 2D, similar to that observed in ARPES experi- 
ments on copper oxides. 



II. POSSIBLE TYPES OF GREEN'S FUNCTION 
RENORMALIZATION. 

Let us start with some qualitative discussion of pos- 
sible manifestations of NFL behavior. Green's function 
of interacting system of electrons is expressed via Dyson 
equation (in Matsubara representation, e n = (2n+l)irT, 
£ P = v F (p-p F )) as 1 : 
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(i) 



In the following, we shall use rather unusual definition 
of renormalization ("residue") Z - factor, introducing it 



Q 

via . 



or 



G(e n ,£, p ) — Z(e n ,£ i p)GQ(e n ,£,p) — — (2) 



^( £ T1) £p) 



fan -£p)G(s n ,£ p ) 
(3) 

Note that Z{e n ,^ p ) is in general complex and actually 
determines full renormalization of free - electron Green's 
function Go(e n ,^ p ) due to interactions. At the same 



1 Despite our use of Matsubara representation, below we consider 
e n as continuous variable. 
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time, it is in some sense similar to standard residue renor- 
malization factor used in Fermi liquid theory. 

Let us consider possible alternatives for Z(e n ,£ p ) be- 
havior. 



A. Fermi liquid behavior. 

In normal Fermi liquid we can perform the usual ex- 
pansion (close to Fermi level and in obvious notations), 
assuming the absence of any singularities in E(e„,p): 



S(e„,C P ) « E(0,0)+ie t 



<9E(e„,£ p ) 



d{ie n 



<9E(e„,£ p ) 



(4) 

In the absence of static impurity scattering E(0, 0) is real 
and just renormalizes the chemical potential. Then we 
can rewrite Q as: 

1 7 

G(s) = - 



■>■£,. 



<9E 
d{i£n) j o 



as 

^pJo 



where we have introduced the usual renormalized residue 
at the pole: 

<9E 



Z = 



1 



1 - 



<9E 



d(is n ) 



d{ie r , 



(6) 



and spectrum of quasiparicles: 



_9E\ 
d $P Jo 



(7) 



The usual analytic continuation to real frequencies gives 
now the standard expressions of normal Fermi liquid 



theor 



.10.11 



with real < Z < 1, conserving the quasi- 



particle pole of the Green's function. 

In the special case of £ p = 0, i.e. at the Fermi surface 
which is not renormalized by interactions (according to 
Landau hypothesis and Luttinger theorem), we have: 



G(£m £p) 



z 



(8) 



i.e. Z just coincides with the limit of Z(e n — > 0,£ p = 0) 
as defined by J2J), ©, and we have the usual pole, as 
e„ — > 0. Similarly, for e n = 0, we have Z(e n = 0,£ p — > 
0) - Z. 

In general this behavior is conserved not only for the 
case of E(e„,£ p ) possessing regular expansion at small 
e n and £ p , but also for E(e„, £ p ) ~ Max(e™,£p) with any 
a > 1. 



B. Impure Fermi liquid. 

In case of small concentration of random static impuri- 
ties we have E(e„ — > 0, £ p — > 0) — > const, with i?eE(0,0) 



giving again the shift of the chemical potential, while 
imE(0,0) ~ 7, where 7 is impurity scattering rate. For 
the Green's function we have: 



G(£n,£p) 



(9) 



i^n £,p + *7 | e ™| 

so that renormalization factor defined by iJSJ is given by: 



Z( £ „,£ p ) = Z- 



(10) 



For £ p = we just have: 

Z(e n ,C P = 0) = Z 



ie n + Z7 T 



for |e„|->0 



(11) 



while for |e„| < |£ p | 



Z(e n ^0,Q = Z 



i—signe n 

1 



for £ p 







(12) 



i.e. impurity scattering leads to Z - factor being zero at 
the Fermi surface, just removing the usual Fermi liquid 
pole singularity and producing a finite discontinuity of the 
Green's function at e n = 0. This behavior is due to the 
loss of translational invariance of the Fermi liquid theory 
(momentum conservation) because of impurities. In fact, 
Green's function is obtained after the averaging over 
impurity position, which formally restores translational 
invariance, leading to a kind of (trivial) non - Fermi liq- 
uid (NFL) behavior. Note, that this behavior is observed 
for \e n \, |£ p | <C 7, while in the opposite limit we obviously 
have finite Z(e£ p ) ~ Z. 



Superconductors, Peierls and excitonic 
insulators. 



Consider now the case of s - wave superconductor. 
Normal Gorkov Green's function is given by: 



G(e n , £ p ) 



is 



n ~ sp 



( i£n )2-£2_| A | 2 



(13) 



where A is superconducting gap. The same form normal 
Green's function takes also in excitonic or Peierls insu- 
lator, where A denotes appropriate insulating gap in the 
spectrumai. Then: 



% ( e n, £p) 



Max(el,$) 
IAI 2 



{ien? - (U 2 

(z £ „)2 - e P - iap ~ 

^0 for e n ,^ p ^0 



(14) 
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i.e. we have NFL behavior with pole of the Green's func- 
tion at the Fermi surface replaced by zero, due to Fermi 
surface being "closed" by superconducting (or insulating) 
gap. 

Again, Fermi liquid type behavior with finite Z - factor 
is "restored" for |e n |, |f p | > |A|. 

However, the complete description of superconducting 
(excitonic, Peierls) phase is achieved only after the in- 
troduction also of anomalous Gorkov function. Excita- 
tion spectrum on both sides of the phase transition is 
determined by different Green's functions with different 
topological properties^. 

D. Non — Fermi liquid behavior due to interactions. 

Non - Fermi liquid behavior of Green's function due to 
interactions may appear also in case of singular behavior 
of £(e„,£ p ) — > oo for e n — ► and £ p — > 0, e.g. power 
- like divergence 2 of £(£ n ,£ p ) ~ Max(e~ a , £~ Q ) with 
a > 0. Obviously, in this case we have Z(e n — > 0,£ p — ► 
0) — > 0, and we again have zero of the Green's function 
at the Fermi surface. 

Another possibility is singular behavior of derivatives 
of self - energy in Q, e.g. in case of X(e„,£ p ) ~ 
Max(s", £p) with < a < 1, leading to weaker than the 
usual pole singularity of Green's function at the Fermi 
surface. 

Both types of behavior are realized within Tomonaga - 
Luttinger model in one - dimension 1 —, where asymptotic 
behavior of G(is n , £ p ) in the region of small £ p ~ e n can 
be expressed as: 

G (Zn ~ Cp) ~ ^7 (15) 

with a' < 1/2. For a' > 1/2: 

G(e„~£ P )~ A + Bsl*'- 1 (16) 

For 3/2 > a' > 1: 

G{e n ^^ p )^A + Be n + Ce 2 n a '- 1 , etc. (17) 

with the value of a' determined by the strength of inter- 
action. 

Special case is the so called "marginal" Fermi liquid 
behavior assumeS^ for interpretation of electronic prop- 
erties of CuO-i planes of copper oxides. This is given by: 

£(e n ,fp) ~ Aie„ln l — (18) 

UJ C 

where A is some dimensionless interaction constant, and 
lo c is characteristic cut - off frequency. If we formally use 



2 Additional logarithmic divergence can also be present here! 



(0) at finite e„, we obtain: 

UJ C 

In this case "residue at the pole" of the Green's func- 
tion (Z-factor) 3 goes to zero at the Fermi surface it- 
self, and again quasiparticles are just not defined there 
at all! However, everywhere outside a narrow (logarith- 
mic) region close to the Fermi surface we have more or 
less "usual" quasiparticle contribution — quasiparticles 
(close to the Fermi surface) are just "marginally" defined. 
At present there are no generally accepted microscopic 
models of "marginal" Fermi liquid behavior in two - di- 
mensions. 



III. FLUCTUATING GAP MODEL. 

Physical nature of FGM was extensively discussed in 
the literaturei*2iiMi£i£iL2iii. It is based on the picture of 
an electron propagating in the (static!) Gaussian ran- 
dom field of (pseudogap) fluctuations, leading to scat- 
tering with characteristic momentum transfer from close 
vicinity of some fixed scattering vector Q. These fluctu- 
ations are described by two basic parameters: amplitude 
A and correlation length (of short - range order) 
determining effective width k = of scattering vector 
distribution. 

In one - dimension, the typical choice of scatter- 
ing vector is Q — 2pp (fluctuation region of Peierls 
transition)***, while in two - dimensions we usually 
mean the so called "hot spots" model with Q = 
{■n/a,ir/a^^. These models assume "dielectric" (CDW, 
SDW) nature of pseudogap fluctuations, but essentially 
the same formalism can be used in case of superconduct- 
ing fluctuations^. 

The case of superconducting (s - wave) pseudogap fluc- 
tuations in higher dimensions is described actually by the 
same one - dimensional version of FGM?^S. 

Attractive property of models under discussion is the 
possibility of an exact solution achieved by the complete 
summation of the whole Feynman diagram series in the 
asymptotic limit of large correlation lengths £ — > ooa&. 
In case of finite correlation lengths we can also perform 
summation of all Feynman diagrams for single - elec- 
tron Green's function, using an approximate Ansatz for 
higher order contributions both in one dimension 5 and in 
two dimensions 2,3 . Similar methods of diagram summa- 
tion can be also applied in calculations of two - particles 
Green's functions (vertex parts)2iiMiLiiiii. 



3 Note that 1191 . strictly speaking, can not give correct definition 
of the "residue", as standard expression JSJ is defined only at 
the Fermi surface itself, where 1191 just does not exist. Thus in 
what follows we prefer rather unusual definition given in 01 . 
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Our aim will be to demonstrate that nearly all aspects 
of NFL behavior discussed above can be nicely demon- 
strated within different variants of FGM. 



A. One — dimension. 

We shall limit ourselves here only to the case of in- 
commensurate pseudogap (CDW) fluctuations^. Com- 
mensurate case£*& can be analyzed in a similar way. Note 
that the same expressions apply also for the case of su- 
perconducting (s - wave) fluctuations in all dimensions. 

In the limit of infinite correlation length of pseudogap 
fluctuations we have the following exact solution for a 
single - electron Green's function^ii: 



ie n + £, 




r 2 + £2 



for e n -> 0, Z P - (20) 



where Ei{—x) denotes integral exponential function and 
we used the asymptotic behavior Ei{—x) ~ \n(^'x) for 
x — > (In 7' = 0.577 - Euler constant). Then, using 
we immediately obtain: 




(21) 



Precisely the same result is obtained if we define for finite 

£n,£,p- 



Z(s n , £p) — 



1 



(22) 



d(ie n ) 



similar to JBJ. Note that due to |e n | -C A, |£ p | <C A we 
obviously have Z > 0, but the usual pole of the Green's 
function at the Fermi surface ("point") of the "normal" 
system is transformed here to zero due to pseudogap fluc- 
tuations. Because of topological stability 9 , the singular- 
ity of the Green's function at the Fermi surface is not 
destroyed: zero is also a singularity (with the same topo- 
logical charge) as pole. But actually FGM gives an ex- 
plicit example of a kind of Luttinger or "marginal" Fermi 
liquid with very strong renormalization of singularity at 
the Fermi surface. 

Consider self - energy corresponding to Green's func- 
tions HDJ: 



ISr, 



f 2 



CA 2 



(23) 



so that taking, for brevity, £ p = and e n — > we get: 



£( £ „^0,£p = 0) 



Icr. 



A 2 



1 



is 



n 1 / 1 e n 



30 



(24) 



i.e. the divergence of the type discussed above. 

In case of finite correlation lengths £ = kT 1 of pseudo- 
gap fluctuations we have to use continuous fraction rep- 
resentation of single - electron Green's function derived 
in Ref£ to obtain renormalization factor as (e n > 0): 



+ Sp + >V F K 



+ Hp + 3iv F * 



(25) 



which can be studied numerically. 

In Fig. ^we show typical dependences of renormaliza- 
tion factor Z(e n ,^ p ). In all cases it goes to zero at the 
( "bare" ) Fermi surface and pole of Green's function dis- 
appears. Essentially, this strong renormalization starts 
on the scale of pseudogap width, i.e. for |e„| < A and 
|£p| < A, reflecting non - Fermi liquid behavior due to 
pseudogap fluctuations. 




FIG. 1: Typical dependencies of Z(e n , £ P ) - factor in one 
- dimensional FGM with finite correlation lengths: depen- 
dences of Z(e n — 0, £ p ) and Z(e n ,£,p = 0) on e„ and £ p for 
wfk/A = 0.1. At the insert: dependences of ReZ(e n = 0,£ p ) 
on £ p for different values of k (in units of ^). Both e„ and 
£ p are given in units of A. 



However, the role of finite correlation lengths £ (finite 
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k) is qualitatively similar to static impurity scattering 4 
and more detailed calculation shows, that the behavior 
of Z - factor at small e„ <C vfk, and <C vpK (with 
e n > 0) is as follows: 



^ for £ 



A 



0, 



(26) 



with a(vpK/A) — > for k — > 0, as seen from Fig. In 
terms of Green's function this behavior corresponds to: 



G(e n , £p) 



_L ( Vf k \ £ " j 
A V A ) ie n 



V 



1 (VpK 



i^al — I (27) 



Thus, for finite k, there is no zero of Green's function for 
e n = and £ p = 0, it remains finite as in impure system. 




FIG. 2: Dependence of a 



' VpK \ 

, A J 



on inverse correlation length. 



Vanishing of renormalization factor Z(e n ,£ p ) at the 
"bare" Fermi surface is in correspondence with general 
topological stability arguments^ - in the absence of static 
impurity - like scattering the pole singularity of the 
Green's function is replaced by zero. In the presence 
of this additional scattering this zero is replaced by finite 
discontinuity, i.e. singularity still persists. 



B. "Hot spots" model in two — dimensions. 

In two dimensions we introduce the so called "hot 
spots" model. Consider typical Fermi surface of electrons 
moving in the CuOi plane of copper oxides as shown in 



Fig. If we neglect fine details, the observed (e.g. in 
ARPES experiments) Fermi surface (and also the spec- 
trum of elementary excitations) in CuOi plane, in the 
first approximation are described by the usual tight - 
binding model: 

e(p) = — 2t(cosp x a + cosp y a) — At' cosp x a cosp y a (28) 

where t is the nearest neighbor transfer integral, while t' 
is the transfer integral between second - nearest neigh- 
bors, a is the square lattice constant. 
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FIG. 3: Fermi surface in the Brillouin zone and "hot spots" 
model. Magnetic zone appears e.g. in the presence of anti- 
ferromagnetic long - range order. "Hot spots" correspond to 
intersections of its borders with Fermi surface and are con- 
nected by the scattering vector of the order of Q = — ). 

Phase transition to antiferromagnetic (SDW) state in- 
duces lattice period doubling and leads to the appearance 
of "antiferromagnetic" Brillouin zone in inverse space as 
is also shown in Fig. |3| If the spectrum of carriers is 
given by l|28|l with t' = and we consider the half - 
filled case, Fermi surface becomes just a square coincid- 
ing with the borders of antiferromagnetic zone and we 
have a complete "nesting" — flat parts of the Fermi sur- 
face match each other after the translation by vector of 
antiferromagnetic ordering Q = (±7r/a, ±7r/a). In this 
case and for T = the electronic spectrum is unstable, 
energy gap appears everywhere on the Fermi surface and 
the system becomes insulator, due to the formation of 
antiferromagnetic spin density wave (SDW) 5 . In the case 
of the Fermi surface shown in Fig|31 the appearance of 
antiferromagnetic long - range order, in accordance with 
general rules of the band theory, leads to the appearance 
of discontinuities of isoenergetic surfaces (e.g. Fermi sur- 
face) at crossing points with borders of new (magnetic) 
Brillouin zone due to gap opening at points connected by 
vector Q. 



4 This is due to our approximation of the static nature of pseudo- 
gap fluctuations. 



J Analogous dielectrization is realized also in the case of the for- 
mation of the similar charge density wave (CDW). 
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In the most part of underdoped region of cuprate phase 
diagram antifcrromagnetic long - range order is absent, 
however, a number of experiments support the existence 
of well developed fluctuations of antiferromagnetic short 
- range order which scatter electrons with characteris- 
tic momentum transfer of the order of Q. Similar ef- 
fects may appear due CDW fluctuations. These pseu- 
dogap fluctuations are again considered to be static and 
Gaussian, and characterized by two parameters: ampli- 
tude A and correlation length £ = In this case 
we can obtain rather complete solution for single - elec- 
tron Green's function via summation of all Feynman di- 
agrams of perturbation series, describing scattering by 
these fluctuation o 1 ! 2 !? . This solution is again exact in 
the limit of £ — ► oc£, and apparently very close to an 
exact one in case of finite Generalizations of this ap- 
proach for two - particle properties (vertex - parts) are 
also quite feasible. 

We shall start again with an exact solution for £ — * oo 
(or k — 0) 2 . First, let us introduce (normal) Green's 
function for SDW (CDW) state with long - range order 
(see e.giii): 



(i£n - £p)(«£n - £p-Q) - W 2 



(29) 



where W denotes the amplitude of SDW (CDW) periodic 
potential and £ p = e(p) — \i. Then we can write down 
appropriate Z - factor as: 



Z{e n , £ p ) 



(*en-&)(fen-60- W 2 



(30) 



where we have denoted for brevity: £ P = £i and £ p -q = 
6- In the following we shall be mainly interested in the 
limit of £„ — > and £1 — > 0, i.e. on the approach to 
the "bare" Fermi surface. Note that 6 = defines the 
so called "shadow" Fermi surface. We have 6 = £2 = 
precisely at the "hot spots" . In the following it is conve- 
nient to introduce a complex variable: 

z = (ie n - ii){ie n - 6) (31) 

which becomes small for e n) 6>6 — * 0. 



1. Incommensurate combinatorics. 

In case of incommensurate (CDW) pseudogap fluc- 
tuations, an exact solution for the Green's function of 
FGM in the limit of correlation length £ — > 00 takes the 
form similar to 1201^ and we get (averaging l|30|) with 
Rayleigh distribution for W): 



Z(z) 



00 2W 

dW—^e'^ — 
A 2 z 



W 2 



f°° dC < z z . ( z \ 



Then, for z->0we get: 
Z{z -> 0) 



A 2 



hi 



(33) 



At the "bare" Fermi surface we have £1 = and in the 
following we limit ourselves to e n > 0. Then, from l|33|l 
we can easily find limiting behavior of Z{z). Just quoting 
some results we have: 



1. For £„ < |&|: 



2 A 2 

i.e. "impure" - like linear behavior in e n . 



(34) 



2. For e n |6I (i-e. also at the "hot spot", where 
6 = 0): 

ReZ(e n » |6Ui = 0) « -||m (V||) +^f> (35) 

i.e. (for £2 = 0) NFL behavior similar to one - 
dimensional case. 

Note that we always have ImZ = for £2 = 0, i.e. at 
the "shadow" Fermi surface and in particular at the "hot 
spot" itself. 



2. Spin - fermion combinatorics. 

Consider now spin - fermion (Heisenberg) model for 
pseudogap (SDW) fluctuations 2 . In this case we again 
obtain FGM, but with gap distribution is different (from 
Rayleigh distribution) and instead of (|32|l we have: 



'2tt 



W 2 

dW —e 



(¥) 



V^Jo (f) 3/2 z -( 



w 2 



1 



«c — —e 



r(3/2) (-z) 3 / 2 



2tt (Ai) 3/2 



exp 



2(f) 



Thus, for z — > we obtain: 
2r(3/2) 



Z(z) 



2(f), 



3/2' 



2(f) 



r(-i/2) 



2(f) 



Then on "bare" Fermi surface (£ p = 0) we have: 
Z{e n ^ 0,6,^ = 0) = 



(36) 



2r(3/2) 



-r(-i/2) 



£n(£n + «6) 



'A 1 ' 



+ 



£n(£n + j6J 

2(f) . 



3/2' 



(37) 
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In particular, for £2 = 

Z(e»->0,6=6 
_ T(3/2) ej 

^ (f ) 



we have JmZ = and: 
= 0) = iteZ(e B -»0,6 = 6 = Q) 



(38) 



so that we obtain quadratic NFL behavior of Z - factor. 
Again let us present some results on limiting behavior: 

1. For e„ < |6| : 

i?eZ(e„ « 161,6 =0) = 



2F(3/2) 



2(f) 



£n|6l 



3/2' 



(39) 



i.e. NFL "zero" behavior. 



2. For £„ 3> |6| (i.e. also at the "hot spot" where 
6 = 0): 



ReZ(e n ^> 6,6 = 0) 



r(3/2) 



■ A 2 s 
v 3 / 



(40) 



i.e. again NFL "zero" behavior. 



In the general case of finite correlation lengths £ = re -1 
we have to perform numerical analysis using the recur- 
sion relations proposed in Refsi2*2L. Again we use the basic 
definition of Z - factor given in To calculate self - 
energy E(e n ,£ p ) of an electron moving in the quenched 
random field of (static) Gaussian spin fluctuations with 
dominant scattering momentum transfers from the vicin- 
ity of characteristic vector Q, we use the following recur- 
sion procedure 2 * 3 - which takes into account all Feynman 
diagrams describing the scattering of electrons by this 
random field. The desired self-energy is given by 



S(e n ,£ p ) — £/c = i(£ n ,£ p ) 



(41) 



with 6. = e(p) - fi (cf. (H) and 

E fc (e B ,£ p )=A 5 



s(k) 



ie n + fi- e fe (p) + inv k K - £ fe+ i(e„, £ p ) 

The quantity A again characterizes the energy scale of 
pseudogap fluctuations and k = is the inverse corre- 
lation length of short range SDW fluctuations, efc(p) = 
e(p + Q) and v k = |«p + q| + |^p + Cjl for odd k while 
£ k(p) = e(p) and v k = \Vp\ + \v^\ for even k. The velocity 
projections Vp and are determined by usual momen- 
tum derivatives of the "bare" electronic energy dispersion 
e(p) given by (|28l) . Finally, s(k) represents a combinato- 
rial factor with 



s(k) = k 



(43) 



for the case of commensurate charge (CDW type) fluctu- 
ations with Q = (n/a, ir/a) 5 . For incommensurate CDW 
fluctuations 5 one finds 



s(k) 



fc+i 

i 

2 



for odd k 
for even k. 



(44) 



For spin - fermion model of Ref^, the combinatorics of 
diagrams becomes more complicated. Spin - conserving 
scattering processes obey commensurate combinatorics, 
while spin - flip scattering is described by diagrams of 
incommensurate type ("charged" random field in terms 
of Ref. 2 ). In this model the recursion relation for the 
single-particle Green function is again given by gg, but 
the combinatorial factor s(n) now acquires the following 
form 2 : 



s(k) 



fc+2 



for odd k 
for even k. 



(45) 



Below we only present our results for spin - fermion com- 
binatorics, as in other cases we obtain more or less similar 
behavior of renormalization factors. 

In Fig. 0]we show the results of numerical calculation 
of ReZ(s n , 6< = 0) at different points taken at the "bare" 
Fermi surface, shown at the insert. For comparison we 
show data obtained in the limit of infinite correlation 
length £ — > 00 (or k = - exactly solvable case) and for 
finite na = 0.01 (i.e. £ = 100a). It is clearly seen that in 
both cases ReZ ~ 1 at "nodal" point D, except at very 
small values of e n , while in the vicinity of the "hot spot" 
(points A and C), and also at the "hot spot" itself (point 
-B), ReZ becomes small in rather wide interval of e n < A. 
This corresponds to more or less "Fermi liquid" behavior 
for "nodal" region (vicinity of Brillouin zone diagonal), 
with a kind of "marginal" Fermi liquid or Luttinger liquid 
(NFL) behavior as we move to the vicinity of the "hot 
spot" . 




FIG. 4: Dependence of ReZ on e n (in units of t) at different 
points Fermi surface (corresponding to t = — OAt and fi = 
— 1.3t) in "hot spots" model (spin - fermion combinatorics 
of diagrams) with correlation lengths £ — > 00 (re = 0) and 
£ _1 a = rea = 0.01. Pseudogap amplitude A = O.li. At the 
insert we show the "bare" Fermi surface and points, where 
calculations were done. 

For completeness in Fig. [5]we show similar comparison 
of dependences of ImZ on e n at the same characteristic 
points on the Fermi surface and for the same parame- 
ters as in Figs. 0J It is only important to stress once 
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so that spectral density in the case of long - range 
(CDW,SDW) order has the following form: 

A w {e,^) = --ImG R (e 1 ^) = 

= (e - - fr)(e - &) ~ TF 2 ] Siff n(2 £ - - 6) 

(47) 

Accordingly, for FGM with correlation length (^oowe 
have: 

POO 

A(e,£ p ) = / dWV w A w (e,£ p ) (48) 
Jo 

where 7-V is distribution function of gap fluctuations, 
depending on combinatorics of diagrams and leading to 
the following separate cases, already considered (or men- 
tioned) above: 



FIG. 5: Dependence of ImZ on e„ (in units of transfer in- 
tegral t) at different points Fermi surface (corresponding to 
t' = — 0.4t and fi — — 1.3t) in "hot spots" model with finite 
correlation length £ -1 a = Ka = 0.01 (spin - fermion combi- 
natorics of diagrams). Pseudogap amplitude A = O.lt. At 
the insert — the "bare" Fermi surface and points, where cal- 
culations were done. 



again, that only at the "hot spot" itself (point B) we 
have ImZ — 0, so that Z becomes real, and shows de- 
pendence on e n more or less equivalent to that proposed 
for "marginal" Fermi liquids (or Luttinger liquids). 

In all cases we observe vanishing of renormalization 
factor Z(e n ,Zp) at the "bare" Fermi surface. In the ab- 
sence of static impurity - like scattering due to finite val- 
ues of correlation length £ = kT 1 the pole singularity of 
the Green's function is replaced by zero, reflecting topo- 
logical stability of the "bare" Fermi surface (Luttinger 
theorem)^. In the presence of this scattering, singular- 
ity of the Green's function at topologically stable "bare" 
Fermi surface remains in the form of finite discontinuity. 



C. Spectral density and Fermi surface 
"destruction" in "hot spots" model. 

Let us return to (|29[) and perform the usual analytic 
continuation to real frequencies: ie n — > e + iS. Then we 
obtain: 



1. Incommensurate combinatorics. 

In the case of incommensurate CDW - like pseudogap 
fluctuations we have: 

2W _wl 

V W = -£2 e ^ ( 49 ) 

- Rayleigh distribution^^. Then, from 148|) we obtain: 

£ — £ 2 _ (E-£l)(f-£2) 
= ~^~ e 23 X 

x0[(e-£ 1 )(e-6)]« 5 n(2E-fc-$,) 

(50) 

For £ = 0we have: 

A(e = (U P ) = ^e-^Ofotdsignfa + b) (51) 
For £i — * ±0 we get: 

A(e = 0, £ p -> ±0, 6) = ± J^(±6) (52) 

so that within the Brillouin zone A(e — 0, £ p ) is nonzero 
only in the space between "bare" Fermi surface and 
"shadow" Fermi surface. This qualitative result is con- 
firmed below, for all other combinatorics, for the case of 
"pure" FGM with £~ 1 = K = 0. 

2. Commensurate combinatorics. 



6 





(e - £i)(e - 6) - + iS(2e - & - &) 



(46) 



In the case of commensurate CDW - like pseudogap 
fluctuations we have&: 



V w 



1 



-e 2A^ 



(53) 
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Gaussian distribution. Then, from i|48|) we obtain: 

s-i2 



1 



(^-Sl)(»-g2) 

-e x 



(54) 



with the same qualitative conclusions as in incommensu- 
rate case. 



3. Spin - fermion combinatorics. 



1 

f 


A=0.1t 
,6=0.01 











In the case of SDW - like pseudogap fluctuations of 
(Heisenberg) spin - fermion model^ we have gap distri- 
bution: 



Then, from (|48|l we obtain: 



2\3/2 



(55) 



i n J- \ I 7^~\ (»-tl)(e-;2) 

= 1 V(e-6)(e-6) c x 

V2tt (^) 3/2 
xO[{e - fr)(e - &)]sign(2e - 6 - 6) (56) 

again with the same qualitative conclusions as in incom- 
mensurate case. 

For the general case of finite correlation lengths £ = 
spectral densities can be directly computed using 
analytic continuation of recursion relations (|4 1 1> - (|42H to 
real frequencies^. 

Actually, two - dimensional contour plots of A{e — 
0,£ p ) (which are in direct correspondence with ARPES 
intensity plots) can be used for "practical" definition of 
renormalized Fermi surface and provide a qualitative pic- 
ture of its evolution in FGM with the change model pa- 
rameters 6 . 

In Fig. El we show typical intensity plots of spectral 
density A(e = 0, £ p ) in Brillouin zone for the "hot spots" 
model both for the case of infinite correlation length 
= k = and for finite (large!) correlation length 
£~ 1 a = na = 0.01 (spin - fermion combinatorics of dia- 
grams, in other cases behavior is quite similar) and for 
different values of pseudogap amplitude A. We see that 
these spectral density plots give rather beautiful quali- 
tative picture of the "destruction" of the Fermi surface 
in the vicinity of "hot spots" for small values of A, with 
formation of typical "Fermi arcs" as A grows, which is 
qualitatively resembling typical ARPES data for copper 
oxidesiML. 





FIG. 6: Intensity plots of spectral density A(e — 0, £ p ) in 
Brillouin zone for the "hot spots" model (t' = — OAt and \x = 
— 1.3t) for the case of infinite correlation length = k = 
and for finite correlation length = na = 0.01 (spin - 

fermion combinatorics of diagrams) with different values of 
pseudogap amplitude A. "Bare" Fermi surface is shown by 
dashed line. 



D. Superconducting d — wave fluctuations. 



As we noted above, the case of superconducting s - 
wave pseudogap fluctuations simply reduces to one - di- 
mensional FGM. Much more interesting is the case of 
superconducting d - wave fluctuations in 2D. 

To obtain exact results for the case of infinite cor- 
relation length = k = we have only to make 
simple replacements in the above expressions for the 
"hot spots" model with incommensurate combinatorics: 
£2 — * — £,1 = —£,p and A — * A p , where A p defines the 
amplitude of fluctuations with d - wave symmetry: 



A p = —A [cos^a) — cos(p y a)] 



(57) 



Note that for the free electrons A(e = 0, £ p ) = <5(§ p ), so that 
appropriate intensity plot directly reproduces the "bare" Fermi 
surface. 



with A now characterizing the energy scale of pseudogap 
fluctuations. 

Then (|31fl reduces to z — — (e^ + £p) and for Z-factor 



10 



we immediately obtain an expression, similar to (|21|l : 

ei+e P ( ei+e P \ „.( 4+g \ 

,-f^infyf^Uo 
a p V a p / 

for e n -» 0, i v -> (58) 

again replacing the pole singularity by zero at the "bare" 
Fermi surface, except the "nodal" at the diagonal of the 
Brillouin zone, where A p = (cf. 

Instead of l|50(l . we get spectral density as: 

s 2 -i 2 

A (e, C P ) = £ -^e~^f6{E 2 e v )signs (59) 

which is nonzero only for |£ p | < e. As a result, for e = 
we have A(e — 0,£, p ) = for A p ^ 0, and it is different 
from zero only at the intersection of Brillouin zone diag- 
onal with "bare" Fermi surface, where A p given by (|57|) 
is zero. At Fermi surface itself we have: 

A(e,t p =0) = ^e- j k (60) 

with two maxima at e = ±A/V2. 

Considering the general case of finite correlation 
lengths £ = we again perform numerical analysis 
using the recursion relations introduced for this problem 
in Ref. 3 , using the basic definition of Z - factor given in 
©. To calculate self - energy £(£„,£ p ) of an electron 
scattered by static fluctuations of superconducting order 
parameter with d - wave symmetry, we use the following 
relation (similar to H42|) h sligthly generalizing relations 
derived in of Ref. 3 : 

Sfc(e„,£ P ) = 

= A pS (fc) 

ie n - (-l) fc £ P + ikn{\v*\ + K|) - E fc+1 (e n ,£ p ) 

(61) 

where s(k) is defined in 144fl . 

In Fig. [7| we show the results for ReZ[e ni £ P = 0) 
again taken at different points of the "bare" Fermi sur- 
face, shown at the insert. Correlation length is £ = 100a 
(na = 0.01) and A = O.li. It is clearly seen that ReZ = 1 
precisely at the "nodal" point D (where A p = 0), but in 
other point on the "bare" Fermi surface ReZ is strongly 
renormalized in rather wide intervals of e n < |A p |, going 
to zero with e n — > 0. Thus we again obtain a kind of 
"marginal" Fermi liquid or Luttinger liquid (NFL), but 
qualitatively different from the case of "hot spots" model. 

In Fig. |S| we also show typical intensity plots of spec- 
tral density A(e = 0, £ p ) in Brillouin zone for the case of 
superconducting (d - wave) pseudogap fluctuations with 
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FIG. 7: Dependence of ReZ on e n (in units of t) at different 
points Fermi surface (corresponding to t' — — OAt and fi = 
— 1.3t) in the model of superconducting (d - wave) pseudpgap 
fluctuations with correlation length £ -1 a = ko = 0.01. Pseu- 
dogap amplitude A = O.lt. At the insert we show the "bare" 
Fermi surface and points, where calculations were done. 




Pi Pi 

FIG. 8: Intensity plots of spectral density A(e — 0, £ p ) in 
Brillouin zone (t 1 — — OAt and /i = — 1.3t) for the case of 
superconducting (d - wave) pseudogap fluctuations. Correla- 
tion length £ -1 a = kcl = 0.1 (spin - fermion combinatorics 
of diagrams) for two different values of pseudogap amplitude 
A = 0.3t and A = t. 



correlation length £ _1 a = na = 0.1 and two different val- 
ues of A. We see that these spectral density plots give 
quite different picture of the "destruction" of the Fermi 
surface in comparison with the case of "hot spots" model, 
which also, in our opinion, differs significantly from most 
results of ARPES measurements on copper oxides. Fermi 
surface is sharply defined only in one point (at the diago- 
nal of the Brillouin zone) , where A p given by (|57|l is pre- 
cisely zero, and there are no sharply defined Fermi arcs 
formed close to this point. We observe only some more or 
less wide "dragon - fly wings" formed around this point. 
Note also the absence of any signs of "shadow" Fermi 
surface. 
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IV. CONCLUSION 

We analyzed rather unusual (NFL) behavior of fluctu- 
ating gap model (FGM) of pseudogap behavior in both 
ID and 2D. We studied in detail quasiparticle renor- 
malization [Z - factor) of the single - electron Green's 
function, demonstrating a kind of "marginal" Fermi liq- 
uid or Luttinger liquid behavior (i.e. the absence of well 
- defined quasiparticles close to the Fermi surface) and 
also the topological stability of the "bare" Fermi surface 
(Luttinger theorem). This reflects strong renormaliza- 
tion effects leading to the replacement of the usual pole 
singularity of the Green's function in Fermi liquid by zero, 
thus effectively replacing at the Fermi surface of poles by 
Luttinger surface of zeroes^. In the presence of static 
impurity - like scattering due to the effects of finite cor- 
relation lengths of pseudogap fluctuations this singularity 
is replaced by finite discontinuty. 

In 2D case we discussed effective picture of Fermi sur- 
face "destruction" both in "hot spots" model of dielectric 



(AFM, CDW) pseudogap fluctuations, as well as for qual- 
itatively different case of superconducting d - wave fluc- 
tuations, reflecting NFL spectral density behavior and 
similar to that observed in ARPES experiments on cop- 
per oxides. Intensity plots obtained for the case of AFM 
(CDW) fluctuations are, in our opinion, more resembling 
ARPES intensity data obtained in experiments on copper 
oxides. Note, that this effective picture was also directly 
generalized for the case of strongly correlated metals or 
doped Mott insulators^ using the so called DMFT+'Ef. 
approach of RefiiS. 
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